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We investigate the relative particle number squeezing produced in the excited states of 
a weakly interacting condensate at zero temperature by stimulated light scattering using a 
pair of lasers. We shall show that a modest number of relative number squeezed particles 
can be achieved when atoms with momentum k, produced in pairs through collisions in 
the condensate, are scattered out by their interaction with the lasers. This squeezing is 
optimal when the momentum k is larger than the inverse healing length, k > ko. This 
modest number of relative number squeezed particles has the potential to be amplified in 
four-wave-mixing experiments. 



A macroscopic source of strongly entangled 
atoms would be immensely valuable for such appli- 
cations as Heisenberg limited clocks [1], quantum 
information processing (for reviews see [2]), and 
tests of quantum state reduction [3] . Various meth- 
ods of producing macroscopic entangled atoms us- 
ing Bose-Einstein condensates (BECs) have been 
proposed, such as four-wave-mixing [4] and uti- 
lizing spin-exchanging collisions in a trap [5,6]. 
The goal of this letter is to investigate the rela- 
tive number squeezing produced in excited states 
of a condensate using stimulated light scattering 
[7] . We shall briefly discuss how this relative num- 
ber squeezing can be converted into entanglement. 
The main result of this letter is that if atoms, pro- 
duced in pairs through collisions in the conden- 
sate, are scattered out by their interaction with 
the lasers, then a modest number of almost op- 
timally relative number squeezed particles can be 
produced (as shown in Figure 1). 

We consider a homogeneous condensate of 
atoms whose collisional interaction is described 
by a contact interparticle pseudopotential ^{r) — 
UoS^^^{r). Relative number squeezing in the situ- 
ation considered in this paper originates from pair 
producing collisions in a dilute BEG. To under- 
stand the nature of these pair-producing collisions 
present in a dilute uniform BEG, we first consider 
how weak repulsive interactions affect the ground 
state. Gorrelations between free particle states 
with momenta k and —k (all momenta arc under- 
stood to be vectors) are present in the condensate 
ground state, \(t>o)- These arise from the momen- 
tum conserving interaction between two conden- 
sate atoms colliding that admix momentum states 
k and —k into the ground state. Observation of 
the suppression of the static structure factor at 
low k [9] is evidence for the existence of these cor- 
related pair excitations. In general this homoge- 
neous model should be a good approximation for 
the large condensates as produced in recent exper- 



iments. 

In this paper, we investigate the possibility of 
making these correlations more accessible by trans- 
ferring one particle of such a pair into a state of 
higher momentum, e.g. k into k + Ak, using stim- 
ulated light scattering. We find that squeezing of 
the population difference between the modes —k 
and k + Ak in the final state may be achieved in 
principle only for momenta k on the order of or 
larger than the inverse healing length [10]. 

In the homogeneous model condensate the eigen- 
modes of the weakly interacting Hamiltonian, 
which we define below, are travelling plane waves. 
Furthermore, we impose periodic boundary condi- 
tions in a finite volume V, which makes the mo- 
mentum spectrum discrete. This should be a rea- 
sonable approximation since one can fix the photon 
momenta by tuning the two lasers and the atomic 
momenta by assuming a measurement or an ex- 
traction mechanism to be concentrated on the rel- 
evant momenta. 

We can describe the condensate with a second 
quantized Hamiltonian of the form: 

k fc,/ci,/C2 

(1) 

where Uq — ATrh^a/m, a is the 2-body scattering 
length that characterizes the collisions between the 
atoms, and m is the atomic mass. 

By assuming that there is a macroscopic popula- 
tion in the fc = momentum state and neglecting 
terms of order N^^'^, where N is the total number 
of bosons, we can approximate this Hamiltonian 
to describe a dilute condensate at low temperature 
[11] obtaining 

Hbec ~ 2^ ( + -r^ I «afc + alk<^-k) 
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Note that we have redefined the ground state en- 
ergy to bo zero. Wc then diagonahzc Hbec by 
introducing a hnear transformation, = UkCXk — 
^feOf^fe, from particle operators, a^, to quasipar- 
ticle operators, d/c, where Uk = (1 — 

vk = Pkil - Pir^'\ h = l + yl- Vk{2 + ylfl\ 

Vk = 1^1/1^0 1) ^0 = ii-KaNQ/V, and Nq is the num- 
ber of particles in the fc = momentum state. We 
can write 



H. 
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where Uk = {Tikl/2m)yk\/2 + yj. [10]. Through- 
out this paper we are assuming that our system 
is at, or modestly perturbed from, a quasiparticle 
vacuum or T = 0. The ground state \4>q), for the 
uniform weakly interacting BEC, is then defined 
by dfe|(;io) = for ^ 0, implying 
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This state has the form of a set of two-mode 
squeezed states [12] and consequently has the prop- 
erty of being perfectly relative number squeezed, 
or having zero variance in the relative number of 
particles in k and —k, i.e. [A(nfe — n-fe)]^ = 
{{hk - fi-kY) - {fik - fi-k)"^ = 0. 

Squeezing in the relative number of particles be- 
tween states is not equivalent to entanglement, but 
these squeezed states, if not already entangled, can 
be converted into entangled states. For example, 
a twin Fock state, which is an optimally relative 
particle number squeezed system, can be converted 
into an entangled system via a 50/50 beam-splitter 
[13-15]. The ground state of a dilute uniform BEC 
(eq. (4)) is optimally relative number squeezed. 
However, this is presumably of practical use only 
if the populations of these squeezed states can be 
separated from the condensate ground state. 

We shall quantify the degree of relative num- 
ber squeezing by defining the squeezing parameter, 
ik,k', to be 



Cfe.fe' 
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where nk = {fik) = {o-Xo-k) is the number of par- 
ticles in a momentum state k. This squeezing pa- 
rameter is equivalent to the one introduced by [16] 
in the spin-squeezing formalism to take into ac- 
count the correlations of particles discussed by [17]. 
If the two relevant states are independent and co- 
herent then = 1- On the other hand, for 



optimally relative number squeezed populations of 
the momentum states k and k' , [A(nfc — nfe/)]^ = 
implying ^k,k' = 0. So as (,k,k' decreases, the rela- 
tive number squeezing between the two states im- 
proves. 

Relative particle number squeezed systems can 
be a valuable source for entanglement if the parti- 
cles in the squeezed states can be extracted from 
the condensate. In this letter, we investigate 
how to accomplish this extraction using stimulated 
light scattering. In this method a condensate is il- 
luminated using two laser beams with a wave vec- 
tor difference of AA; and a frequency difference of 
i5 [18]. Since we are approximating quasiparticles 
by travelling plane waves, the outcome of a loca- 
tion measurement of particles is completely unde- 
termined, and the scattering strictly speaking does 
not lead to an extraction, where an atom moves 
out along a defined trajectory. However, since we 
are interested only in populations, we assume that 
particles with a momentum large enough to leave 
a trapped condensate may be described approxi- 
mately by plane waves. In a more realistic for- 
mulation the particles would be described by wave 
packets whose momentum distribution is smeared 
out around the momentum of this plane wave. 

In the remaining part of this letter, we study 
ways to obtain the best possible amount of squeez- 
ing between atomic states of different momenta. 
We calculate the time evolution of the system and 
the relevant populations using perturbation the- 
ory. We treat the effective coupling between the 
laser and the condensate, O, as a perturbation to 
the Hamiltonian (3) describing a weakly interact- 
ing condensate. The results presented below have 
been calculated to the order Vt^ . Hence, we are im- 
plicitly assuming that the lasers do not drive the 
system far out of equilibrium, so that the Bogoliu- 
bov approximation still holds. 

The effective Hamiltonian can be written as 
H = Hbec + Hl, where 
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(subscripts k± are shorthand for k ± Afc). Us- 
ing the quasi-particle transformation already in- 
troduced, the Bogliubov approximation implying 
ao,al — > y/No, and noting that (a^^. -|- a_Afe) = 
{uAk - VAk){oiAk + Q'-Afe) we obtain 
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where Ukk' — UkUk> + ffeUfc', Vkk' = UkVk' + VkUk'- 
Note the factor of (uAfc — ^^Afe) originates because of 
the destructive interference between two resonant 
paths connecting the k = momentum state to the 
momentum states ±Ak [9]. We have approximated 
the internal structure of the atom to be represented 
by a two-level system. The coupling is given by 
n = |r2p/2A, where fl is Rabi frequency of the 
one-photon transition between the internal atomic 
states. The detuning, A, of the laser from the in- 
ternal atomic transition frequency is assumed to be 
sufficiently large to adiabatically eliminate the in- 
termediate state. We have neglected the AC-Stark 
shifts. 



We investigate the squeezing behavior of parti- 
cle populations of modes k ^ ±Ak that are not 
directly coupled to the condensate mode {k = 0) 
by their interaction with the lasers. We find 
that the best squeezing occurs between popula- 
tions of momentum states k + Ak and — fc as 
the two laser beams with their frequency differ- 
ence tuned to S = LJk + u^k+Ak excite a parti- 
cle from a correlated pair in the ground state 
[19]. In this case the interaction Hamiltoninan is 
Hl « -^h,ftvkk+ctl_^_^^&^_^+h..c., which has the 
form of a two-mode squeezing operator. Therefore, 
the gwasiparticle populations at momenta k + Ak 
and —k are optimally relative number squeezed to 
the order of the calculation. Since this simplified 
Hamiltonian is proportional to Vkk+ : it is clear that 
the collisions during the stimulated photon scatter- 
ing are the essential process for the production of 
squeezing. We project these optimally squeezed 
quasiparticle states onto a particle basis to ob- 
tain the particle squeezing parameter £,k+Ak,-k as 
a function of |A:| (assuming Afc = fc/2) for succes- 
sive times as seen in Fig. 1. 




t = 0.01 ms 



t = 0.1 ms 



FIG. 1. The squeezing parameter ^k+Ak,-k, for 
Afc = fc/2, as a function of |fc| for successive times 
for typical experimental parameters: A^o ~ W"^ 
^^Na atoms, Volume V = 10"'^ cm^, a = 2.8 nm, 
Rabi frequency Q = (27r) 1.8 MHz, and detun- 
ing A = (27r)lGHz implying no = 10"cm"^, 
nko/2m = (27r) 1.5 kHz and n/(?ifco/2m) = 1.0. Note 
that the squeezing improves with increasing time and 
increasing momentum. 

As time increases the squeezing improves 
and the squeezing parameter asymptotically ap- 
proaches a minimum value 

^k+Ak-k « + ■vl^)ul^ + {ul + vl)ul 

-2ulul^]/[ul^+ul]. (8) 

Although perturbation breaks down at large times, 
the limit of the validity of our perturbative calcu- 
lation is very close to the asymptotic limit. On 

resonance, the particle populations of atoms with 
momenta k + Ak and —k are approximately given 

by 

Uk+Ak ~ ^^fc+ + i {^tful^ vlk+ > (9) 
n_k^vl + \{Qtfulvlk^, (10) 

respectively. We can see from Figure 1 that the 
relevant states arc approximately relative number 
squeezed {^k+Ak,-k close to zero) in the particle 
regime > ko. For example, if we assume the 
same experimental parameters as used in Figure 1, 
the two lasers' combined pulse with a duration of 
10 ms should produce roughly 60 particles in both 
momentum states 3fco and — 2fco- These momen- 
tum states should be almost optimally squeezed 
with ^fe+Afe,-Afe « 0. 




|A k/kj 

FIG. 2. The squeezing parameter ^Ak,-Ak as a func- 
tion of |Afct for successive times for the same exper- 
imcutal parameters as used in Figure 1. Note that 
tlic squeezing deteriorates with increasing time and 
the squeezing is concentrated at the low momentum 
(|Afc| < fco) states. 



It is also instructive to investigate the lead- 
ing order process in which particles in the con- 
densate mode fc = are directly scattered into 
modes ±Ak. We consider the relative particle 
number squeezing between the momentum states 
k = Ak and k = —Ak. Assuming the laser 
beams' frequency difference is tuned to the res- 
onance Sj= WAfe, the interaction Hamiltonian is 
H]^ ss ^Ufl y/l^{uAk — VAk) ct]\i.+ ii-c., SO there are 
no correlations between quasiparticles to the order 
of our calculation and therefore no relative num- 
ber squeezing between quasiparticlc populations. 
However, when we project the excited quasiparti- 
clc states onto a particle basis there is significant 
relative number squeezing in the phonon (or quasi- 
particlc) regime, \k\ < ko, [20] (as shown in Fig- 
ure 2). This squeezing originates from the initial 
state (cq. (4)) and it is clear that the squeezing 
deteriorates when additional atoms are scattered 
from the condensate into the modes ±Afc. This is 
contrary to the previous situation where relative 
particle number squeezing originates from the col- 
lisions during the scattering process and improves 
with time. 

In summary, we have considered a situation 
where one can excite relative number squeezed par- 
ticles using stimulated light scattering. We have 
shown that by considering states that are not di- 
rectly coupled to the macroscopically populated 
fc = state, one can in principle excite and extract 
populations of strongly relative number squeezed 
particles. We expect that the dominant process 
to limit this relative number squeezing will be the 
rescattering of atoms in collisions before they leave 
the condensate. Assuming a cubic volume V, we 
estimate the percentage of rcscattered particles to 
be smaller than r = anoV^/^ =42%. We expect 
that the inclusion of such processes into our calcu- 
lation would yield an upper bound for the; achiev- 
able squeezing and an approximate time when this 
optimum is reached. Another measure of the loss 
is the mean decay time due to quasiparticlc col- 
lisions. At zero temperature the main loss mech- 
anism is given by Beliaev damping processes, in 
which cither a quasiparticlc collides with a con- 
densate atom to produce a pair of quasiparticles 
or two quasiparticles collide and thereby produce 
an atom in the condensate and another quasiparti- 
clc. For I A; I ^ ko, the inverse Beliaev decay width 
Tfe = = m/{8Tra^nohk) [21] which may also 
be obtained from kinetic theory, is for our above 
example T2kQ — 3.5 ms, giving an estimate for the 
time when deteriorating processes become signifi- 
cant. 

The particles produced in the relative number 
squeezed modes might be used as seeds for four- 



wave-mixing experiments. Since the introduction 
of uncorrelated particles can be minimized, the 
squeezing in the output of such experiments can 
be greatly improved. Note that a similar ampli- 
fication method has recently been demonstrated 
[22] for the production polarization-entangled pho- 
tons, using stimulated emission inside a nonlinear 
medium. 
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